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What can one do with a given tunable quantum device? We provide complete symmetry criteria 
deciding whether some effective target interaction(s) can be simulated by a set of given interactions. 
Symmetries lead to a better understanding of simulation and permit a reasoning beyond the limi¬ 
tations of the usual explicit Lie closure. Conserved quantities induced by symmetries pave the way 
to a resource theory for simulability. On a general level, one can now decide equality for any pair 
of compact Lie algebras just given by their generators without determining the algebras explicitly. 
Several physical examples are illustrated, including entanglement invariants, the relation to unitary 
gate membership problems, as well as the central-spin model. 
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I. INTRODUCTION 


Thanks to impressive progress on the experimental 
side, many small- and medium-scale quantum devices are 
now ready for applications ranging from quantum metrol¬ 
ogy M to quantum simulation With quantum 

information processing as one of the driving but long¬ 
term goals (e.g., UM3), one of the pressing questions 
is, what can one do with these devices now? This prob¬ 
lem clearly falls into the remit of quantum systems and 
control engineering , an area naturally receiving increased 


interest |13H15| both experimentally and theoretically. 

Control theory offers a well-known characterization 
of the operations that a quantum device is capable of 
on Lie-algebraic grounds [lil. fl6M2lj | . In this work, we 
simplify the question to the Hamiltonian membership 
problem of (finite-dimensional) quantum simulation. It 
amounts to deciding, for a set of given control interac¬ 
tions V , whether a set of effective target interactions Q 
can be simulated—without having to establish controlla¬ 
bility via nested (and hence tedious) commutator calcula¬ 
tions for the so-called Lie closure. Our results reduce the 
Hamiltonian membership problem to the straightforward 
solution of homogeneous linear equations. 

In the setting of the controlled Schrodinger equa¬ 
tion [22| (taken as a bilinear control system 17, 


2[/(f) = [-iHt + Y^ v=2 ~ ™ v (t)H v \ U(t), (1) 

we ask whether the given set V := {iH \,..., iH p } of in¬ 
teractions (which may include a drift term) generates an 
effective interaction iH p+ i or more generally any inter¬ 
action from a set Q := {iH p+ i,..., iH q } assuming all 
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H v are represented by Hermitian matrices henceforth. If 
so, then for every evolution time r > 0 of a simulated 
interaction iHk 6 Q, there is a solution U (t) of the sim¬ 
ulating system © for 0 < t < 8 and controls u v (t) such 
that V generates a unitary U(6) = exp(— irHk) in the 
simulation time 8 starting from the identity at t = 0 
0 H3T13 ■ M ■ I n this sense, Hamiltonian simulation of 
a particular Hamiltonian Hj. can be considered as an in¬ 
finitesimal version of creating a particular unitary gate. 
It also generalizes the universality (or full controllabil¬ 
ity) question of whether all Hamiltonians can be simu¬ 
lated (or equivalently whether all unitary gates can be 
obtained) [3, [28} H2T1I1- In the context of gates, a fa¬ 
miliar elementary example is that all unitary gates in 
an n-qubit system can be obtained [32} by combining lo¬ 
cal gates with CNOT gates. However, the approach of 
the pioneering age of decomposing every target gate into 
a sequence of CNOT and local gates is, in practice, all 
too often imprecise or slow. So implementing gates or 
simulating Hamiltonians with high fidelity rather asks 
for optimal control techniques, as explained in a recent 
roadmap [42|. As a precondition, here we step back to 
the Hamiltonian level and give criteria for simulability 
and controllability. 


II. MAIN IDEA 

We solve the decision problems of simulability (and 
controllability) by just analyzing the symmetries of the 
Hamiltonians of given setups. We show that this decision 
requires considering both linear and quadratic symme¬ 
tries, where linear symmetries of a Hamiltonian H com¬ 
mute with H, while quadratic symmetries of H are those 
commuting with the tensor square ( iH <g> 1 + 1 ® iH). 
The term quadratic symmetry is motivated, since the 
tensor square generates XJ <8> U just as iH generates the 
unitary U. 

More precisely, our goal is to get a symmetry-based 
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understanding of how a set V of available interactions 
can simulate a set Q of desired effective quantum inter¬ 
actions in the sense that the Lie closures coincide, i.e. 
(V) = {VL> Q). We circumvent brute-force calculation of 
the Lie closure not only because high-order commutators 
can entail a significant growth in the appearing matrix 
entries and may lead to instabilities in numerical com¬ 
putations, but first and foremost because it provides no 
deeper insight into the problem. Our symmetry anal¬ 
ysis leads to a much more systematic understanding of 
Hamiltonian simulation and quantum system dynamics 
in general. It provides a powerful argument to decide 
under which conditions a desired Hamiltonian can, in 
fact, be simulated, or in turn, which explicit simulations 
or computations are impossible in a given experimental 
setup. 

Let us summarize our line of thought: As short-hand, 
let the linear symmetries of V (analogously for any set 
of matrices) be expressed via the commutant V' which 
consists of all matrices S £ C dxd that commute (i.e ., 
[5, iH v \ = 0) with each element iH u £ C dxd 0 f V f43t |. 
Obviously, for Q to be simulable by V, it is necessary 
that Q may not break but rather has to inherit the 
symmetries of V, so Ami\P'\ = dim[('PuQ)']. However, 
a complete symmetry characterization is nontrivial. It 
rather requires the following two steps: The first is to 
introduce quadratic symmetries [l9j as those linear sym¬ 
metries of the system artificially doubled by the tensor 
square V® 2 := {iH v ®t d +t d ®iH u for v £ {l,...,p}}. 
It defines the quadratic symmetries by its commutant 
•p( 2 ) ._ y_ Secondly, let C denote the center Q 
of the commutant (VLiQ)' and consider the central pro¬ 
jections of V and VUQ onto C. With these stipulations, 
we summarize our main result: 

Main result (see Result {7\ below): The given interac¬ 
tions V simulate the desired interactions Q in the sense 
(V) = (VUQ) if and only if V and Q share the same 
quadratic symmetries (i.e., dimfT 7 ^ 2 ^] = dim[('PUQ)^ 2 ' ) ]) 
[condition (A)] and the central projections of V and VUQ 
onto C are of the same rank [condition (B)]. 

Let us emphasize that our approach goes beyond the 
ubiquitous use of linear symmetries in physics, since lin¬ 
ear symmetries provide only an incomplete picture of 
Hamiltonian simulation. The application of higher sym¬ 
metries is the key here. It is interesting to note that 
essentially only the quadratic symmetries (and no higher 
ones) in condition (A) are necessary to characterize the 
dynamics of a quantum system. One obtains a complete 
description together with the auxiliary condition (B). 

Some remarks also summarizing known approaches are 
in order. The quadratic symmetries are stronger than the 
linear ones; actually they include them and thus condi¬ 
tion (A) implies that the linear symmetries also agree. 
Example [T] below illustrates why matching the linear 
symmetries does not suffice to ensure simulability. As 
shown in a companion paper one can decide if a 
subalgebra f) C g of a compact semisimple Lie algebra 
g actually fulfills fi = g (e.g., (V) = (PLiQ)) just by 


analyzing quadratic symmetries. But Example [2] eluci¬ 
dates why condition (A) alone does not, in the general 
compact case, imply simulability. Only after fixing the 
central projections by condition (B) the quadratic sym¬ 
metries decide simulability. 

On a much more general scale, condition (B) closes 
the gap to completely characterizing equality in f) C g 
now for all compact Lie algebras (generated by skew- 
Hermitian interactions) beyond the semisimple ones of 
ilhl . Simplifying within the Lie-algebraic frame, our sym¬ 
metry approach to decide simulability and the member¬ 
ship Q C ( V) can thus be seen as a major step bey ond the 
well-established Lie-algebra rank condition U- BOj 
and beyond the limited first use of quadratic symmetries 
to establish full controllability in fl9| . 


III. SYMMETRIES 

In this section, we elaborate our method and estab¬ 
lish necessary and sufficient conditions for Hamiltonian 
simulation to arrive at Result [T] below. We also describe 
important properties of linear and quadratic symmetries 
and discuss two illustrating examples. Example [T] high¬ 
lights the importance of quadratic symmetries for decid¬ 
ing Hamiltonian simulation and their relevance for en¬ 
tanglement invariants. The necessity for the auxiliary 
condition (B) is made evident in Example [5] 

The linear symmetries of A 4 C C dxd are identified [l9j 
with the commutant AT given as 

AT := {S £ C dxd s.t. [S, M] = 0 for all M £ M}. 

The commutant includes all complex multiples of the 
identity l d and it forms a vector space of dimension 
dim(AT). A smaller set of matrices typically shows more 
symmetries, i.e., for AT C AT one has AT A AT and 
AT = AT iff dim(AT) = dim(AT). By Jacobi’s iden¬ 
tity (i.e., [S,[M 1 ,M 2 ]] = [[M 2 ,S\,M 1 ] + [[S,M 1 ],M 2 ]), 
any symmetry S that commutes with both M\ and M 2 
also commutes with their commutator [Mi, M2]. So, A4 
and the Lie algebra (A4) it generates have the same com¬ 
mutant: AT = AT 2 if (AIi) = (AT). 

In our context, this implies that iH p+ \ cannot be sim¬ 
ulated by V unless V = (V U {iH p+ 1 })', i.e., coinciding 
symmetries are a necessary but not sufficient condition. 
This is because the converse does not hold as the follow¬ 
ing basic example illustrates: 

Example 1: The pair interaction iH zz := iZ\Z 2 cannot 
be simulated by the local interactions V = {iX\,iY^, 
iX 2 , iY 2 } of a two-qubit system [Hj in spite of coinciding 
(trivial) commutants V = (V U {iH zz })' = CI 4 . 

Thus, we further discuss quadratic symmetries [IqJ de¬ 
fined by the commutant to the tensor-square [47) . 

AT 2) := (A! 02 )' = {S £ C d2xd2 such that 

[S, M®l d +l d ®M] = 0 for all M £ M Q C dxd }. 
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FIG. 1. (Color online) Visualization of Example [T] (a) No 
linear symmetries besides the identity exist for both the fully 
controllable system and the local interactions, (b) The dou¬ 
bled system reveals nontrivial quadratic symmetries corre¬ 
sponding to separate permutations (1,3) and (2,4). 

The tensor-square commutant always contains (the sub¬ 
space spanned by) the identity l d 2 and the SWAP or com¬ 
mutation matrix [48| . Also, the quadratic symme¬ 
tries include all linear ones, i.e., 5j ®lrf+lrf<g>5j £ 
for S\ £ M' . And by Jacobi’s identity [5l|, one finds 
(AJi)^ = (M 2 ) (2> if (Adi) = (M 2 ). As above, in our 
context this implies that iH p+ i cannot be simulated by 
P unless P ^ = (P U {iH p+ i})^ holds. 

Example Q] (completion): The relevant tensor-square 
commutants have different dimensions dim[p( 2 )] = 4 

and dim[(P U {iH zz })^] = 2, so iH zz cannot be simu¬ 
lated. Naturally, (V U {iH z z })^ contains li 6 and the 
commutation matrix K^, which is related to the joint 
permutation (1,3)(2,4) of tensor components in C 16x16 , 
while V contains two additional quadratic symmetries 
related to the separate permutations (1,3) and (2,4); 
see Fig. [Q Evidently, the local interactions of P can¬ 
not generate entanglement. Hence, a quadratic symme¬ 
try in P^ 2 * has a physical interpretation as an entangle¬ 
ment invariant. Indeed, the concurrence p32j of a two- 
qubit pure state | ip) can be defined as [(V’l(V’|li 6 — M(i, 3 ) 
—M(2 i 4)+M( 1j3 )(2 i 4)|'!/’)|'!/’)] 1 ^ 2 /2 HM 1 , where the ma¬ 
trix M p is defined by the permutation p. Any quadratic 
symmetry S £ V relates to a degree-two polynomial 
invariant Tr[p®pS l ] in the entries of the density matrix p 

m- 

Remarkably, symmetries beyond quadratic ones (i.e., 
those of the tensor square) are not required for a neces¬ 
sary and sufficient condition for simulability [58|. Con¬ 
cerning the tensor-square commutant, we build on two 
important classification-free results of [45j for compact 
Lie algebras |H, [ 6 (j (as generated by skew-Hermitian 
matrices iH v \. For (PUQ) being semisimple (and com¬ 
pact), Ref. [45| first shows that (P) = (PUQ) holds if and 
only if dimfP^ 2 )] = dim[(PuQ)*' 2 '*]. Beyond the semisim¬ 
ple case, any compact Lie algebra g can be uniquely de¬ 
composed as g = s ® c into its semisimple part s and 
its center c (where s := [g,g] and [g, c] = 0 [i3|). So 
Ref. [45j secondly verifies that the semisimple parts of (P) 
and (PUQ) have to agree if dim[p( 2 )] = dim[(PUQ) (2) ]. 
When generalizing from semisimple to arbitrary compact 
Lie algebras, the equality of the two tensor-square com¬ 
mutants implies that (P) and (PUQ) agree-except for 
the central elements (commuting with all the other ones). 
These commuting interactions require condition (B) to 
fix the central projection thus resulting in the following 


complete characterization: 

Result 1: Consider two sets V := {iHi,... ,iH p } and 
Q := {iH p+ 1 ,..., iH q } of (skew-Hermitian) interactions, 
and let C a denote elements of a linear basis spanning 
the center C of the commutant (PUQ)L For the central 
projections, define the matrix T by its entries T a p := 
Tt(C) y iHp) for 1 < a < dim(C) and 1 < /3 < q as 
well as T by T af 3 := Tr (C^iHp) for 1 < (3 < p. Then 
V simulates Q in the sense (V) = (PUQ), if and only 
if both conditions (A) dim[p( 2 )] = dim[(7 : ’UQ)^ 2 ' 1 ] and 
(B) rank(T) = rank(T) are fulfilled. 

Condition (B) of Result [T| is a basic linear-algebra test 
solely based on linear symmetries. Each of the matrices 
T and T depends on both V and Q. In Example [Q iH zz 
could not be generated as condition (A) is not satisfied. 
Before proving Result [T] the following example provides 
a helpful illustration of condition (B): 

Example 2: In a two-qubit system, consider a dipole 
coupling combined with a tilted magnetic field, i.e., V := 
{i(2Z 1 Z 2 -X 1 X 2 -Y 1 Y 2 ), i(X 1 -Y 1 +X 2 -Y 2 )}. We inves¬ 
tigate whether a Heisenberg-type interaction of the form 
Q a := {i(X 1 X 2 +Y 1 Y 2 +Z 1 Z 2 )j or one particular inter¬ 
action of pairing type (i.e., Qb := {i(X 1 Z 2 +Z 1 X 2 +Y L Z 2 
+ZiY 2 )}) can be simulated. Condition (A) is satisfied in 
both cases as the quadratic symmetries of V, VU Q a , and 
VUQb all coincide (there are 16 of them). The three lin¬ 
ear symmetries also agree. Moreover, with the mutually 
commuting operators 
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/10 0 o\ 

[ 0010 „„ r i 
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forming a basis of the commutants V = ('PU Q a )' = 
(PU Qb)', they also span the (three-dimensional) center C. 
For the central projections, one thus gets the matrices 


/ 0 0 0 \ / 0 0 0 \ ~ ~ 

T “ = (n-r)’ Tb = (^n)’ andr “ =Tb 



Condition (B) reveals rank(T a ) ^ rank(T a ), rank(T&) = 
rank(T{,), so Q a cannot be simulated by V, while Qb can. 
Note the isomorphy types of (V), (PU Q a ), and (PUQh) 
are su(2)®u(l), su(2)®u(l)®u(l), and su(2)®u(l). 

Proof of Result Q] Decompose the compact Lie alge¬ 
bras (P) and (PUQ) into their semisimple parts and cen¬ 
ters. If condition (A) holds, the semisimple parts coin¬ 
cide, (P) = s + c-p and (PUQ) = s + c-pyQ with tv Q 
c-pu Q- Take the unique decomposition iH 1 = iHf + iH\ 
with iHf £ s and iHf £ c-puQ- Since [iH n ,iH m ] £ s, the 
real-linear span of {iHf,... ,iHf} agrees with c-p, while 
the one of {iHf ,..., iHf} equals Cpuc- It follows that 
c p = cp\jQ iff the dimensions of the two real-linear spans 
agree; this in turn is equivalent to the dimensions of the 
two complex-linear spans being equal because all relevant 
Lie algebras are compact (see Corollary 1 of Theorem 1 
in Chapter IX, Sec. 3.3 of ®). 

The center c of a compact Lie algebra g lies within 
the center C of its matrix commutant g': [c,g\ = 0 for 
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c £ c, g £ g implies c C g'; likewise [s,c] = 0 for s £ g' 
shows that c C C. Given the basis { C a } of C, introduce 
its dual basis {B a } with respect to the Hilbert-Schmidt 
scalar product via Ti -(C^Bp) = 8 a ,p- So any C € C can 
be written as C = Y) a Tr(C^C) B a . Define the matrix 
K entrywise by K a p := Trfor 1 < a < dim(C) 

and 1 < /3 < q, and similarly K by K a p := Tr (C^iHp) 
for 1 < (3 < p. Hence the dimension of cp u q agrees with 
the rank of K , and the dimension of c-p equals the rank 
of K. 

Now, for any b £ s, there are two elements bi,b 2 £ 5 
with b = [ 61 , 62 ] (Chapter I, Sec. 6.4, Proposition 5 
of (HI). Th us Tr(C a 6 ) = Tr(C a [b 1 ,b 2 ]) = Tr(C a 6 A) 

- Tr^aMi) = Tr(C Q M 2 ) - Tr(b 2 C a bi) = Tr(C' a 6 1 6 2 ) 

— Tr(C a 6i6 2 ) = 0, where the third equality follows as 

C a £ C and b 2 commute, and cyclic permutations in the 
trace imply the fourth equality. Moreover, Tr(CT6) = 
-TrtCi&t) = —Tr (C a b) = 0. Hence^ Tr (CliHp) = 
Tv{C} y iH c p) 1 which implies T = K and T = K. In sum¬ 
mary, c-p = c-puQ iff the ranks of T and T agree, which 
proves Result [0 □ 


IV. ALGORITHMICS AND BEYOND 

Both linear and quadratic symmetries can readily 
be computed by standard linear algebra: Linear sym¬ 
metries S £ C dxd are determined by the commutant 
and can be obtained by solving the linear equations 
(ld<8>M—M t <g>ld) vec(S) = 0 jointly for all M £ M C 
<C dxd [H [50]. Here, vec (S) is a column vector of 
length d 2 stacking all columns of S [4^|. The dimen¬ 
sion of the solution is d 2 —r where r denotes the rank 
of the matrix formed by vertically stacking the matri¬ 
ces td®M—M t ®td- Likewise, the quadratic symmetries 
S £ <C d xd (given by the tensor-square commutant) 
just amount to solving [l d 2 (gi(Af(g)l d +l (i (g)M) — 
J rtd®M) t ®t d p] vec(S) = 0 jointly for all M £ M. The 
preceding discussion explains how to explicitly determine 
linear and quadratic symmetries. This allows us to test 
condition (A) (i.e., dim^ 2 )] = dim[('PUQ)^]) by com¬ 
paring the dimensions of the quadratic symmetries for V 
and VUQ. 

As the commutant (PUQ)' represents the linear sym¬ 
metries of PL)Q, its center C is readily obtained by solving 
the linear equations (td®M— vec(C) = 0 and 
(ld®S—5* ( 8 ) Id) vec(C) = 0 jointly with M extending 
over all M £ PL)Q and S over all S £ (VL)Q)'. Solv¬ 
ing for C yields a basis C a of the center C, and one can 
determine the matrices T and T as T a p = Tr 
for 1 < a < dim(C) and 1 < /3 < q, as well as 
T a p = Tr (C^iHp) for 1 < /3 < p. Since condition (B) 
is given by rank(T) = rank(T), it can easily be tested by 
elementary linear-algebra computations comparing the 
ranks of T and T. To sum up, Result |T| reduces the 
Hamiltonian membership problem to straightforward so¬ 


TABLE I. Central-spin model of Example [3] number n of 
spins, Lie dimensions dim^T 3 }) = dim((7 3 UQ)), the isomor- 
phy type, dimensions of quadratic and linear symmetries (i.e. 
dimpP (2) ] = dim[(7 3 UQ) (2) ] and dim[P'] = dim[(PUQ)']), and 
ranks of the central projections [i.e., rank(T) = rank(T)]. 


n 

Lie- Isomorphy No. of symmetries 

dim. type quad. lin. 

Rank of 
c. proj. 

case (a): Jk ~ 1 



2 

15 su(4) 

2 1 

0 

3 

38 su(2)©su(6) 

8 2 

0 

4 

78 su(4)©su(8) 

50 5 

0 

5 

137 su(2)©su(6)©su(10) 

392 14 

0 

6 

221 su(4)ffisu(8)ffisu(12) 

3528 42 

0 

case (b): Jk — 2 for even k and Jk = 

1 otherwise 


2 

15 su(4) 

2 1 

0 

3 

63 su(8) 

2 1 

0 

4 

158 su(4)ffisu(12) 

8 2 

0 

5 

396 su(2)ffisu(6)©su(6)©su(18) 

32 4 

0 

6 

796 su(4)©su(8)fflsu(12)©su(24) 

200 10 

0 


lutions of homogeneous linear equations. 

Example 3 (central-spin model): Consider a central 
spin interacting with n— 1 surrounding spins via a star¬ 
shaped coupling graph (where the surrounding spins may 
be taken as uncontrolled spin bath) ItiTl foil . The inter¬ 
actions amount to a drift term (tunneling plus coupling) 
and just a local Z-control on the central spin, V := {iX\ 
+i Jk(XiXk+YiYk+ZiZk), iZi}. We ask whether 

the central spin can be fully controlled, i.e., if Q := {iXi} 
can be simulated. Depending on the interaction strengths 
A £ M for k > 2, different cases are possible: (a) with 
Jk = 1 and (b) with Jk = 2 for even fc, and Jk = 1 
otherwise. 

Computational results for the central spin model have 
been obtained using exact arithmetic [641 ] for a moder¬ 
ate number of spins as detailed in Table |T] These results 
vary significantly for different coupling strengths Jk . But 
our approach for deciding simulability allows for analytic 
reasoning even beyond specific choices of Jk- For Hamil¬ 
tonian simulation, it thus provides a powerful technique 
to analyze and understand the dynamics of general quan¬ 
tum systems. This even holds if the symmetries cannot 
be calculated explicitly. Showcases for the strength of 
explicit symmetries are given in Examples [l][3j while Ex¬ 
ample [3] also makes use of symmetries implicitly (in parts 
where they cannot be calculated explicitly) via the proofs 
of the Appendix [Aj These proofs motivate the following: 

Conjecture: In the central-spin model of Example [3j 
the central spin is fully controllable for a finite number 
of spins and any choice of Jk (i.e., iX\ can be simulated, 
and the surrounding spins can be uncoupled by applying 
the control). 
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V. DISCUSSION 

Similar to the Hamiltonian membership problem for 
interactions solved here, one may address membership 
for groups, e.g., (i) in the (prototypical) discrete case, 
(ii) in connected compact Lie groups, and (iii) in non- 
connected compact groups including finite groups. 

In discrete groups (i), asking the question (a) if Q = 
{U p + 1 } is (exactly) contained in the group generated by 
the unitaries V = {f/i,..., U p } is undecidable for SU(iV) 
(at least for N > 4) 65| . Yet the question (b) of approx¬ 
imate universality [66l l67l|. i.e., if all unitaries in SU(iV) 
can be approximated, is decidable [65|, [68j by compar¬ 
ing the matrix algebra generated by elements U„®U V for 
with its equivalent for SU(IV) (plus other condi¬ 
tions). Still, the tedious algebra closure is needed, similar 
to the Lie closure. Question (b) is equivalent to compar¬ 
ing the topological closure of the group generated by V 
to SXJ(N) and thus leads to (ii). 

In continuous groups (ii), Result |T] applies to decide if 
two connected, compact Lie groups (given by their in¬ 
finitesimal generators) are equal: 

Result 2: Given two sets V and Q of (skew-Hermitian) 
interactions, the elements of V simulate the ones of Q 
and vice versa iff both (V) = (VUQ) and (Q) = (PUQ) 
hold, where each condition can be tested by Result [T| 

Our findings do not generalize to non-connected com¬ 
pact groups (iii), nor are they implied by the represen¬ 
tation theory of compact groups. In particular, finite 
groups with trivial quadratic symmetries [S', U„(&U V \ = 0 
only (known as group designs [69[) do not contradict our 
work. 


VI. CONCLUSION 

We have presented a complete symmetry approach to 
decide Hamiltonian simulability, i.e., whether given drift 
and control Hamiltonians can simulate a target (effective) 
Hamiltonian in finite dimensions. Quadratic symmetries 
lead to an understanding that allows one to algebraically 
prove simulability in classes of many-body systems where 
the usual computational assessment via the Lie closure 
is infeasible. This is exemplified by proving simulability 
for interesting cases of the central-spin model (see the 
Appendix 0 for which only very restricted cases were 
addressed before [63j. 

Achievability of specific target interactions is particu¬ 
larly important for fault-tolerance, where the simulation 
of a particular Hamiltonian (or universality) is needed 
only on logical subspaces and not globally. While linear 
symmetries have often been used in those cases [70l - [73j . 
going a step further by applying quadratic symmetries to 
ensure controllability or simulability on a noise-protected 
subspace could be an interesting application, simplify¬ 
ing com plic ated system-algebraic analysis. For instance, 
in Ref. [741 ] . we examined standard scenarios of noise- 


protected subspaces, where controllability was (moder¬ 
ately) easy to assess. However, in more realistic settings, 
analyzing quadratic symmetries and their restrictions to 
protected subspaces is anticipated to be much easier than 
establishing Lie closures over restricted subspaces. 

Moreover, our results on quadratic symmetries distin¬ 
guishing local properties from global ones can be gener¬ 
alized into an overarching framework that encapsulates 
concurrence (Example [1]) and links naturally to entangle¬ 
ment detection via a quadratic invariant of the quantum 
system under local transformations in [7bl - l78i |. 

Our findings imply that for any nonsimulable interac¬ 
tion, a related resource is lacking. In Example 1 it simply 
was entanglement, but more generally we can character¬ 
ize lacking resources as induced by conserved quantites 
arising from quadratic symmetries. This paves the way 
toward a resource theory of quantum simulability. 
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Appendix A: Central-spin controllability for 
different levels of generality 

In this Appendix, we analyze under which conditions 
on the coupling coefficients Jk the central spin is con¬ 
trollable in Example [3] We collect proofs for this con¬ 
trollability under varying assumptions. Recall the set 
V = {iH\, HI 2 } of control interactions, where iHi = iX 1 
+ * 2 Jk(XiX k +YiY k +ZiZ k ) and iH 2 = iZ\, as well 

as the target interaction Q = {iX{\. Assuming that con¬ 
dition (A) of Result [T| holds, there exists an element iC in 
the center of PC{iX{\) such that iXi+iC £ (V). Since 
iZi £ (' P ), one obtains — [iZi, [iZ\, iXi+iC}]/A = iX L £ 
(V). Thus, it follows that it suffices to verify condition 
(A) in the different cases below, which is equivalent to 

,4 

showing that D(iXi)v = 0 holds for all vectors v £ C 
with d := 2 n and D{iHi)v = D{iH 2 )v = 0. Here, the 
linear operator 

D(M) ■.=[t d 2®(M®l d +l d ®M) 

~(M®l d +l d ®M) t ®l d 2\ £ C d4xdi 

is a shortcut in order to define the linear equations 
D{M)v = 0 for the matrix M £ C dxd and quadratic 
symmetries S where v := vec(S). One naturally 
obtains that both of the equations [D(Mi), D(M 2 )] 
= D([Mi, M 2 ]) and exp[D(Mi)]D(M 2 ) exp[— D(Mi)] = 
D[exp(Mi)M 2 exp(—Mi)] hold for all matrices Mi, M 2 . 
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Proposition 1: The interaction iX\ can be simulated 
if all couplings J k are either (a) equal, i.e., J k = J , (b) 
equal up to an odd integer o k , i.e., J k = Jok where Ok 
may depend on fc, or (c) Q-linear independent. 

Proof. Consider the definitions iH zz := i ]C” =2 J k Z\ Zk 
= iHi + [ii? 2 ,[iF 2 ,iFi]]/4 6 (V) and iH := iX 1 + 

1 EL 2 JkiXiXk+YiYk) = iHi — iH zz £ (V). We assume 
in the following that D{iH\)v = D(iH 2 )v = D{iH zz )v = 
D(iH)v = 0 holds in order to prove D(iXi)v = 0. 

The joint eigenbasis of the operators D(iZ k / 2) and 
D(iZ\Zk/2) for k £ {2is given by the compu¬ 
tational basis, and its basis vectors are w(b) = \bi) <g> 
••• <8 \b in ) with b k S {0,1}, |0) := (0,1)*, and |1) := 
(1, 0)*. This implies that the eigenvalue equations are 
D{iZ k /2)w{b) = ip k (b)w(b) and D(iZiZ k /2)w(b) = 
iX k (b)w(b), and the corresponding eigenvalues are given 

by 

Hk(b) = i( — Sfc — Sn+/c+S2ra+fc+S3n+fc), 

Xk(b) = 2 ( — SlSfc — S„+lS ra -|_/ c -|-S2n+lS2n+fe+S3n+lS3n+fc)j 

where /z fc (6) £ {-2, -1,0,1, 2}, A k {b) £ {-2, -1,0,1, 2}, 
and Sj := 2bj—l. By checking all of the 2 8 cases 
for s en+ i,s en+ k £ {-1,+1} and e £ {0,1,2,3}, one 
concludes that p k (b) mod 2 = A k (b) mod 2 holds if 
D{iZ\)w(b) = 0. Recall that D{iZ\)v = D(iH zz )v = 
0 and expand v as v = 'Yf b abw(b). It follows that 
the equations D(iZi)w(b) = 0 and D(iH zz )w{b) = 0 
hold for ab ^ 0 as each w(b) is an eigenvector of 
D{iZi) and D{iH zz ) = J2 k =2 2 Jk D(iZ 1 Z k /2). Assum¬ 
ing D(iZi)w(b) = 0, this also means that the relation 
p z (b) mod 2 = A zz (b) mod 2 holds for the eigenvalue 
ifi z (b) of D{iZk/2) and the eigenvalue iX zz (b) of 

D{iZ\Zk/2). Moreover, we obtain for ab ^ 0 that 
0 = D(iH zz )w(b) = [iJ2k=2^^kXk{b)]w{b) and, conse¬ 
quently, 2 2.J k X k {b) = 0. 

The proof depends now on the particular cases, and we 
prove in each case that fi z (b) mod 2 = 0: For the case (a) 
with J k = J, it follows that A zz (b) = Y ^ k =2 ^k(b) = 0. 
This implies that fi z (b) mod 2 = 0. In case (b), we 
obtain J k = Jo k and A zz (b) mod 2 = X^l -=2 W(&) mod 

2 = J2 k =2 °kX k {b) mod 2 = 0, which also shows that 

p z (b) mod 2 = 0. For case (c), Jf,k =2 2 </fcAfc(&) = 0 means 
that A k (b) = 0 for all k since the couplings J k are Q- 
linear independent and A k (b) £ Z. In particular, it fol¬ 
lows that A zz (b) = Afc(6) = 0, which proves again 

that p z {b) mod 2 = 0. 

Define the operator 

n 

W := expbr^T D(iZ k /2)\. 

k- 2 

Using the properties of J2k=2 D(iZ k /2), one gets that 
the equation Ww(b) = w(b) = w(b) holds for 

each w(b) with ab ^ 0, where the last equality fol¬ 
lows from p z (b) mod 2 = 0. Thus, we obtain Wv = 


WJ2 b abw(b) = J2b a bWw(b) = J2 b a b w {b) = v. We 
also have that WiD^Xk+Y^W^ = ■iD[G{X- [ X k 
+yiF fc )Gt], using the notation 

n n n 

G := exp(7r^ iZ k /2) = exp(niZ k /2) = iZ k . 

k =2 fc=2 k=2 

It follows that WiD(XiX k +YiY k )W^ = [}}{', fe „ =2 

(iZk^Xk+Y^kX-iZk,)} = ^iD{X l Xk+Y 1 Y k ) since 
Z k X k Z k = -X k and Z k Y k Z k = -Y k . Naturally, 
WD{iX\)W^ = D(iX i) is also satisfied. 

One can now verify that 

0 = WD(iH)v = WD(iH)W*Wv 

n 

= [DiiXi) -Y,i J k D ( x i x k+Yi Y k)] v 

k—2 

n 

= D(iH)v - 2 ^2 iJ k D(X 1 X k +Y 1 Y k )v. 

k= 2 

This implies ’^f J r ] i=2 iJkD(XiXk+Y{Yk)v = 0, and one 
concludes that D(iH)v — J2 k -2 i JkD(Xi X k +Yi Y k )v = 
D(iX i)« = 0. ^ □ 

The techniques in the proof of Proposition [T] can be 
generalized in order to establish the following result: 

Proposition 2: The interaction iX\ can be simulated if 
J k = J for 2 < k < no and J k = 2 J for no < k < n. 

Proof. We establish again all the properties of the first 
two paragraphs in the proof of Proposition [T] Then, it 
follows that YJkL 2 J^k{b) + Y2=n 0 +i 2 ^fe( & ) = Let 
i/J^z\b) be the eigenvalue of J2 k =2 D{iZ k /2). One ob¬ 
tains that /il °\b) mod 2 = 0 for each w(b) with a b ^ 0. 
Define the operator 

n 0 

W (0) := exp[7r^ D{iZ k / 2)]. 

k~2 

We apply the properties of Y^kL 2 D(iZ k /2) and con¬ 
clude that the equation = e 1 ^w(b) = 

w(b) holds for each element w(b) satisfying a b ^ 0, 
where the last equality follows from p> z ^ (b) mod 2 = 
0. Thus, we obtain W^v = J2 b a b w (b) = 

J2b otbW^w{b) = J2b a b w (b) = v. We also have that 
W^iD{X 1 X k +Y 1 Y k ){W^ = iD\GW{X x X k +YiY k ) 
(g(°))t] U si n g the notation 

no no no 

G (0) := exp(7r^ iZ k /2 ) = exp(-K iZ k /2) = iZ k . 

k—2 k—2 k—2 

It follows that 

WWiD^Xk+YiYkXW^)' 

n 0 

= iD[ (iZ k ,){X 1 X k +Y 1 Y k )(-iZ k „)] 

k',k"=2 
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= -iD^Xk+Y^k) 

if 2 < /c < 77-0 since Z^X^Z^ — X^ nnd Zj^Yf^Z]^ — 
and W<+.D(X 1 X fc +Y 1 Y fc )(VF( 0 + = iD(X 1 X k +Y 1 Y h ) if 
k > no- Naturally, D(iXi)(W^)^ = D(iX i) is also 
satisfied. 

One can now verify that 

0 = W w D{iH)v = W^ 0) D{iH){W {0) )^W {0) v 

n 0 

= [D{iX{) - ^+JT>(X 1 X fe +y 1 y fc ) 

k =2 
n 

+ ^ 2iJD(X 1 X k +Y 1 Y k )]v 

k=no~\-l 

n 0 

= D{iH)v - 2 Y iJD(X 1 X k +Y 1 Y k )v , 

k—2 

which implies ]yL ~2 iJfcT^XiXfc+yiYfc)!; = 0. Thus, one 
can conclude that D{iH)v — Y^k =2 iD{XiX k +Y{Y k )v = 
D(iH^)v = 0, where we introduced the notation 

iH^ := %X\ +*Efc =no +12J(x 1 x fe +y 1 y fc ). 

Furthermore, the equation D(iH zz )v = D(iH^)v 
= D(iH 2 )v = 0 implies the important commutator iden¬ 
tity [[[D(iHW),D{iH 2 )],D(iHW)\,D(iH zz )]v = 0. In 
addition, we have 

[[[D(iH^),D(iH 2 )],D(iH^)\,D(iH zz )\ = 

n 

D([[[m i \H 2 ],m%H zz ])=64J 2 3 4 5 6 D(i J2 Y k)- 

k=no~\~l 

Thus, D(i X)fc =no+ i Y k )v = 0. Now, we also get that 

n n 

0=[[D(iH zz ),D(i J2 Y k)],D(i J2 Yk )> 

k=no~\-l k—no~\-l 

= -4 D(iH^)v, 

where iHH 1 := i’52 k =n 0 +l^JZiZk- Considering the ex¬ 
pansion v = '^2b abW (b)i we obtain from = 0 


that the condition A+ (&) = 0 holds for the eigen¬ 
value iXi 1 ] (b) of Y^k= no +i D(iZ\Z k /2) with respect to 
a vector w(b) with a.b ^ 0. As we have for any w(b) 
with D(iZi)w{b ) = 0 that the eigenvalue */+(&) of 
!Cfe=n 0 +i -0(*^fc/2) satisfies A + (&) mod 2 = /+(&) mod 

2, thus we can conclude that /i+(6) mod 2 = 0 for b with 
^ 0. We can now define the operator 

n 

:= exp[7r ^ D(iZ k / 2)]. 

k—TiQ +1 

Using the properties of Efc=n , 0 +1 -D(*^k/2), one gets that 

the equation W^w(b) = e 1 ^w(b) = w(b) holds 
for each w(b) with ctb ^ 0, where the last equality fol¬ 
lows from i^z\b) mod 2 = 0. Thus, we obtain W^v = 
W< ' 1 ' 1 J2b a bUl(b) = J2b a b W< ' 1 ' >W ( b ) = J2b a b w ( b ) = 
v. We also have that W(+.D(XiXfc-|-yiYfc)(W’( 1 + = 
iD[G < ' 1 \XiX k +YiY k )(G ( ' 1 ^ 1 )^] using the notation 

n n 

G (1) := exp(7r E ^^fc/2) — j j 6Xp(7T Z^/q/2) 

A:—no +1 k—TiQ +1 

n 

= p 

/c=n 0 + l 

With these preparations, one can now verify that 
0 = W (1) .D(bff (1) )u = W (1) D{iH)(W (1) yW {1) v 

n 

= [D(iX 1 )-2 J2 iJD(XiX k +YiY k )\v 

k=no~\~l 

n 

= D{iH)v- 4 ^ iJD{X 1 X k +Y 1 Y k )v, 

k-no~\-l 

which hence implies Efc=n 0 +i iJ k D(XiX k +YiY k )v = 0. 
Consequently, one can finally conclude that D{iH^)v — 
ELno+i 2 iJD^Xk+Y^v = 7+X+ = 0. □ 
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